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Abstract 

We consider the hadronic description of the — ► tt + tt~ decay, with the aim 
to investigate the strong phases generated by the final state interactions. The 
derivation of the dispersion relations using the Lehmann-Symanzik-Zimmermann 
formalism and the Goldberger-Treiman method to include inelastic effects in the 
spectral function are presented. We discuss the problem of quark-hadron duality 
and estimate in the hadronic formalism the corrections to the factorized amplitude 
in the heavy quark limit. 
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1 Introduction 



The inclusion of the strong interaction effects in the theory of exclusive nonleptonic B 
decays is a very difficult task. The problem has been investigated recently by many au- 
thors, in particular, for charmless decays into light pseudoscalar mesons, since the strong 
phases of these amplitudes are crucial for the determination of CP-violating phases in 
present and future experiments JTJ. The first measurements of the branching fractions of 
the B decays into tttc and ttK final states |@]-[f§], considerably stimulated the theoreti- 
cal and phenomenological work devoted to these processes in various approaches. In the 
so-called "naive factorization approximation" ||, the matrix elements of the operators en- 
tering the weak effective hamiltonian are expressed as products of meson decay constants 
and hadronic form factors, which are evaluated in a phenomenological way. An obvious 
defficiency of this approximation is the renormalization scale dependence of the results, 
expressed as /i-dependent Wilson coefficients multiplied by //-independent hadronic form 
factors. Improvements to the factorization approximation were discussed in several papers 
Recent calculations of the B — > tttt decay amplitude were performed either in the 
generalized QCD factorization approach P|-||14|, or by more conventional perturbative 
QCD methods pl-pf. 



The nonleptonic B decays were also investigated recently in a hadronic approach, in 
which a part of the strong dynamics accompanying the weak decay is described using the 
unitarity of the S'-matrix, dispersion relations and Regge phenomenology JL7|-[^7|. In the 
present paper, we apply this approach to the particular case of B® — > 7r + 7r~ decay. One 
aim of our study is to compare the predictions of the hadronic and the partonic treatments 
and to test the validity of quark-hadron duality. 

In the next section, we discuss the derivation of dispersion relations with respect 
to the momentum squared of the external particles, by applying the Lehmann-Symanzik- 
Zimmermann (LSZ) reduction formalism |28j to the S'-matrix element of the decay process. 
In Sect. |3| we explain the Goldberger-Treiman procedure to include inelastic contributions 
in the spectral function and apply it to the amplitudes of the decay B® — > n + ir~ . In Sect. | 
we consider the problem of quark-hadron duality, and estimate in the hadronic formalism 
the corrections to the factorized amplitude produced by the final state interactions in the 
heavy quark limit. Section |5] contains our conclusions. 



2 Dispersion relations for the decay amplitude 



We consider the decay amplitude 

A(B° d - 7T+7T-) = (Tr+ih) ir-(k2),out\H w (0)\B° d (p),m) , (1) 

where the "in" and "out" states are defined with respect to the strong interactions and 
7i w is the weak effective hamiltonian density 
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7= E 

V j=u,c 
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In this relation, 0« are local AB = 1, AS* = operators, and Cj the corresponding Wil- 
son coefficients, which take into account perturbatively the strong dynamics at distances 
shorter than Using the expression (0) of the weak hamiltonian, the decay amplitude 
(0) can be splitted in two terms 

A(B° d - 7T+7T-) = V ud V* b A u + V cd V; b A c , (3) 

with the CP violating phase 7 = Arg(V^* fe ) appearing in the first term. 

The physical amplitude ([I]) is calculated for p = k\ + k 2 at on-shell values of the 
momenta, p 2 = m 2 B , k\ = m 2 , k\ = m 2 . The extrapolation to off-shell external momenta 
can be achieved by the LSZ reduction formalism P8| . In Refs. we applied this 

technique to the expression ([J) of the amplitude. As we shall see below, it is more 
convenient to start from the S-matrix element 

S B o^ + „- = (-K + (kx)7r-(k 2 ),ont\B d (p),m) , (4) 

where the transition from the "in" to the "out" states includes both the strong and weak 
interactions. The expression (|I]) of the decay amplitude is obtained by expanding the S 
matrix to first order in the weak hamiltonian. However, we can apply the LSZ reduction 
to the B meson in Eq. (f|), which leads to the alternative expression 

A(B° d - 7T+7T-) = —L^+^O 7r-(A: 2 ),out|^o(0)|0) , (5) 
V 2 Po 

where rj B o(x) = JC x (j) B o(x) is the source of the meson B d and (fi B o its interpolating field 
(JC X is the Klein-Gordon operator). We recall that in a Lagrangian theory the source, 
which includes both the strong and weak interactions, can be written formally as 

5(f) B o ^5d^(j) B o 

The matrix element in Eq. @ can be defined for arbitrary s = p 2 = (A;! +A; 2 ) 2 , the physical 
amplitude corresponding to s = m B . We notice that Eq. @ is similar to the definition of 
the pion electromagnetic form factor, where i] B is replaced by the electromagnetic current 
Jn. We can apply therefore the standard methods used in deriving the dispersion relations 



for the pion form factor [29|-[31|. More precisely, by the LSZ reduction of one final meson 



(say, 7r + ) in Eq. (p3), we obtain 



x(7r-(fc 2 )|[^ + (x),r ?B o(0)]|0) - ^g-^ J dxe^5(x ) 
x(Tr~(k 2 )\ikio[(t>^+(x),r] B o(0)] - [d <f) n +{x), r/ B o(0)]|0) , 

where r] w +(x) is the source of the reduced pion. 

The second integral in Eq. ([71) contains equal-time commutators produced by the 
action of the Klein-Gordon operator K, x upon the function 9{xq). As shown in the 
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most general form of this term, called "degenerate", is a constant or a polynomial of the 
Lorentz invariant variables. To calculate the degenerate term, one needs the expression 
of the source t]b in terms of the interpolating field 0^+ or its time derivative, which in 
a hadronic Lagrangian theory might be obtained from the formal expression (0). The 
commutators can be then evaluated in principle by applying the canonical commutation 
rules, satisfied, up to a normalization constant, by the interpolation fields However, 
in the standard model, the hadronic fields are defined in terms of the underlying quark 
and gluon degrees of freedom, and the definition of the off-shell fields might introduce 
ambiguities in the evaluation of the degenerate term (the result depends also on which 
pion, 7r + or 7r~, is reduced, since their quark content is different). 

We now turn to the first term of Eq. (|7]) , which is usually called " dispersive term" , 
and has a more complicated structure as a function of the squared external momenta. 
The integral defines a function holomorphic at the values of the momenta for which it 
is convergent. First, due to the presence of the 9(xq) function, the integral upon xq in 
Eq. (^) is convergent for lmk w > 0, i.e. in the upper half of the k w complex plane. A 
detailed analysis must exploit also the causality properties of the commutator |30| , which 



restricts the integral upon the spatial variables to |x| < \xq\. The difficult part of the 
conventional proofs of the dispersion relations is to show that the integrals upon xq and 
x are convergent for complex values of the external momenta. As discussed in p0 |, it is 



sometimes useful to go to a particular Lorentz frame and consider a particular variable, 
for instance k 10 , instead of trying to think in terms of Lorentz invariants. Also, it is useful 
to treat simultaneously the matrix elements (irn, out\r] B \0), (mi, m\r] B \0) and (7i\t]b\tt), 
which are represented by the same analytic function in various parts of the complex plane 
of the dispersive variable. 

In the present case, it is convenient to choose the system with the unreduced pion 
7r~ at rest (k 2 = 0), when k w = (s — k\ — m 2 )/2m 7T . In what follows we shall either 
work with s variable keeping k\ = m 2 fixed, or with k\ variable at fixed s = m? B . By 
expressing in terms of k\o and the fixed Lorenz invariant momentum, the first term of 
Eq. (|7[) depends only on the variable k W) and is analytic in the complex kio plane, except 
a possible discontinuity along the real axis, given by 

^ ^ 10 2y/4:k 01 po 

x *£5(h + k2-p n )(^~(h)\rh+(0)\n)(n\VBo(0)\0) 

- 5>(A;i +p„)(7r-(A: 2 )|^o(0)|n)(n|77 7r+ (0)|0) 

This expression is obtained formally from (|7]) by replacing i9{xq) by 1/2, inserting a 
complete set of intermediate states in the commutator and using translational invariance 
|30| . In order to evaluate the spectral function, we recall that the sources contain both 



the strong and the weak interactions, the last ones being treated to first order. As we 
shall see below, the spectral function takes different forms, depending on the external 
momentum adopted as dispersive variable. 

Let us assume first that k\ is fixed at the physical value k\ = m 2 , and treat the 
amplitude as a function of the variable s = (ki + k 2 ) 2 . In the reference system chosen 



3 



above kio = (s — 2m 2 T )/2m 7T and = kf — m 2 . It is easy to see that in this case the 
second sum in the expression (B) brings no contribution. Indeed, since k\ = m 2 , the only 



state which contributes is the one-pion state \n) = and (vr^lO) = [pO . 

As concerns the first sum of Eq. (||), the intermediate states n which contribute are 
of two kinds: the first ones are generated by the weak part of the source i]b in the second 
matrix element, and undergo a strong transition to the final state 7r + 7r~, mediated by the 
strong part of rj^. According to Eq. (|3p, the second matrix element is equal to the weak 
decay amplitude of an off-shell meson B, with momentum squared equal to s, while the 
first matrix element is a strong amplitude, evaluated at the cm. momentum squared equal 
to s. Therefore, the contribution of these states in the unitarity sum represents the so- 
called "final state interactions" (FSI). The lowest intermediate state consists of two-pions, 
which produces the lowest branch-point at k w = m n , or, equivalently, s = 4m 2 . 

The intermediate states n of the second type are produced by the strong part of the 
source t]b in the second matrix element, which describes the strong decay of an off-shell 
B meson. The first matrix element, where contributes the weak part of r] n , describes the 
weak transition amplitudes from the state n to the final 7t + tt~ state. These terms are 
usually interpreted as "initial state interactions" (ISI). The lowest physical state which 
can contribute is the pair B*tt, producing the lowest threshold s = (uib* + m^) 2 . 

The whole amplitude can be recovered from the discontinuity by means of a dispersion 
integral. To write it down, we need the asymptotic behaviour of the discontinuity, which is 
difficult to estimate, since it involves off-shell quantities. Assuming, for simplicity, that one 
subtraction is necessary, and combining the possible degenerate terms, discussed above, 
with the subtraction constant of the dispersive part, we express the physical amplitude 
as 

A(B° -> 7T + 7T-) = A(S ) 

m 2 B -s J ^ cr FSI (s) 



7T J - m% - ie){s - s ) 

imi. 



DC 



j d, W! (9) 

(s - m B )(s - sq) 



71 



(m* B +m w ) 2 

where 

Vfsi ~E 5 ( fc i + k 2 -p n )M*(n -> mr)A(B -> n) 

n 

visi » E + fc 2 - Pn)A*{n -> 7i7i)M(B -> n) (10) 

n 

are the spectral functions associated to the final (initial) state interactions, respectively. In 
these relations, A(M) denote the amplitudes of the weak (strong) transitions, respectively, 
evaluated for an off-shell B momentum squared equal to s. Since B is stable with respect 
to the strong interactions, mg < m#» +m 7r , the initial state interactions (the last integral 
in Eq. (^)) do not contribute to the on-shell imaginary part. We mention that a dispersion 



relation similar to Eq. (J^j was derived recently in Ref. |32[ for K — ► tttt decay, starting 
from the definition ([]]) of the decay amplitude, treating the weak hamiltonian 7i w as the 
source of a spurion, and using the Mandelstam representation. 
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As mentioned above, the expression ([?[) can be analytically continued also in the 
variable kf, at fixed s, equal to the physical value s = m\. In this case, in the reference 
system chosen above, fcio = {m 2 B — k\ — m 2 T )/2m n and k 2 = [k\ — (mg + m n ) 2 } [kf — (m# — 
m 7r ) 2 ]/(2m 7r ) 2 . The spectral function is given formally by the same expression (||), but 
now the contributions are different. First, we notice that the second sum in Eq. (|8]), which 
previously vanished on account of k\ = m 2 , now includes terms which are produced by 
the strong and the weak parts of the source rj^, for a variable k\. Just like in the discussion 
above, there are intermediate states which are generated by the strong decay of an off- 
shell pion, and undergo then a weak transition mediated by the weak part of rjB, and also 
intermediate states which are produced by the weak part of r] n and generate afterwards 
a pion and a B meson through a strong interaction. In the first case the lowest branch 
point is at k\ = 9m 2 , corresponding to the intermediate state with three pions, and in 
the second case at k\ = (m# + rn n ) 2 , corresponding to the intermediate state tiB. From 
the connection between kio and k\ in the particular system mentioned above, one can 
see that the branch cut corresponds to negative values of fcio- This means that these 
contributions originate actually from the matrix element (7r jr/^ |7r) , related to the B — > tttt 
decay amplitude by crossing symmetry. 

As concerns the first sum in the spectral function (|8|), the intermediate states which 
bring a nonvanishing contribution have p 2 = (k\ + A^) 2 fixed at the value m\. The strong 
part of t]b gives therefore no contribution, since the lowest state possible (the pair B*ir) 
can not be produced at this energy. On the other hand, the weak part of the source rjs 
brings a nonvanishing contribution, producing intermediate particles which undergo then 
a strong interaction. This contribution represents therefore the final state interactions. 
The delta function implies the condition k\ = (p n — k 2 ) 2 , where p 2 = m 2 B and k\ = m 2 , 
which gives for k\ the allowed range k\ < (ms — m n ) 2 . We notice that, unlike the other 
branch points discussed above, which are determined by the lowest intermediate states in 
the unitarity sum, the range of the variable k\ has a more kinematical nature. 

We express now the whole amplitude in terms of its discontinuity by a dispersion 
integral. Assuming that one subtraction is necessary and including the degenerate terms 
in the subtraction constant, we obtain 

A(B° -> 7T+7T-) = A(k 2 ) (11) 

l m 7T - H) [ J./2 &FSI\k\ ) 



+— — 

7T 








7T 




(ml — 

1 V 7T 


4) 



—oo 
oo 



d^- 



" >1 [kf-ml -ze)(Af 
, 2 vi(K 2 ^ 



9m; 



(k? - m 2 )(k? - k 2 ) 



vr J 1 (k[ 2 - ml){k? - Kg) ' 



where ctfsi is given formally by the same expression (|T0|) given above, evaluated now at 
fixed s = m 2 B and variable k 2 , and 

ax « 5 ( k ± + Pn)A*(n -> irB)M{ir -> n) 

n=3n,... 
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E S(h+p n )M*{n^TTB)A{<K 

n=Bir,... 



n) 



(12) 



We denoted generically, as before, by A(M) the weak (strong) amplitudes, respectively. 

It is easy to see that both dispersion relations (§) and (pTTJ ) lead to the same disconti- 
nuity on shell, equal to the spectral function crpsi evaluated for physical masses. These 
relations might be useful in principle if the decay amplitude can be calculated (by chiral 
theory, lattice, etc) at some particular points (s = s or k\ = Kq), with a better accuracy 
than at the physical points, s = m 2 B and k\ = m 2 , respectively. The complete evaluation 
of the dispersion relations is very difficult, since they involve off-shell quantities in the 
spectral functions. However, as we will show below, the dispersion relation (PD can be 
written in a different form, more convenient for the study of the final state interactions. 
We first remark that the spectral function a FSI defined in Eq. ( ]T0| ) is actually indepen- 
dent of the dispersive variable kf, for fixed (fci + k 2 ) 2 = m 2 B . This can be easily seen by 
performing the phase space integral in the first sum of Eq. (^) in the cm. system, when 
p n = k x + k 2 = and the energy squared is equal to m 2 B . As the weak decay amplitude 
(the matrix element {n\rj B \Qi) / ^/2p G in Eq. (||)) and the invariant strong amplitude (the 
matrix element (tt~ (k 2 ) \rj n + (0) \n) / \/2fcoi) depend both only on the Mandelstam variables 
and the physical masses of the particles involved, the result of the phase space integration 
at fixed s is independent of k\ and contains only on-shell quantities^. 

We notice that the most general form of a function having a branch cut for — fi 2 < 
k\ < (m B — m^) 2 , with a constant discontinuity &fsii is 



.2, n-/2, , a FSI ^ 
7T 



A FS i(k() = V{k[) + 



k\ - {m B - m^) 
H 2 + k 2 



(13) 



where V(kl) is a polynomial (more generally, an entire function), independent on apsi- 
In order to construct the full decay amplitude we must add to the function A FS i the 
contribution of the degenerate terms and the last two dispersion integrals in Eq. (|TT| ) , with 
possible subtractions. All these terms, as well as the polynomial V(kf), are independent of 
the discontinuity ufsi- Therefore, by combining them into a single constant, and choosing 
the scale fi 2 ~ m 2 B , we write the physical amplitude as 



A(B° d - 7T + 7T-) 



A , aFSI 1 

A H In 



7T 



mz 



(14) 



where A is the genuine contribution which remains when the long distance final state 
interactions are switched-off, i.e. ufsi — 0. It is worth mentioning that this separation 
of the final state interactions from the other parts of the dynamics was possible only due 
to the fact that crpsi does not depend on the dispersive variable, allowing us to construct 
Apsi according to (p~3|) . The subtraction constants in the usual subtracted dispersion 
relations, like Eqs. (§) or (|TT|), do not have a similar interpretation: they represent the 
values of the amplitude at some particular points, and depend implicitly on all the spectral 
function in the dispersion relation. 



2 The unusual property of the spectral function to be independent of the dispersive variable fcf , was 
noticed a long time ago 1 30 . 
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We shall end this section with some comments. First, we recall that rigorous analytic 
properties in the external momenta are proved in axiomatic field theory only for a small 
region close to the physical masses [Q. Therefore, the dispersion relations presented 



above can be accepted only as an heuristic conjecture, whose validity remains to be 
tested. We mention also that the off-shell analytic continuation in external momenta is 
in general plagued by ambiguities. They may appear, in the present formalism, in the 
evaluation of the degenerate terms and of the off-shell amplitudes entering the spectral 
functions. Moreover, we notice that even the analytic properties of the off-shell amplitude 
may depend on the specific expression of the on-shell amplitude, used as starting point 
of the extrapolation. For instance, by applying the LSZ procedure to the expression 
(0) of the amplitude, we obtained in Ref. only a part of the dispersive branch-cuts 
written in Eqs. @ and (pTTj) . The contribution of the missing dispersion integrals (namely, 
the FSI contribution in Eq. @ and the ISI contribution in Eq. (|ll])) is hidden in the 
corresponding degenerate terms, which have a different form [pH]. Of course, one expects 
that the amplitude on-shell is recoverd in an univoque way, but the compensation of the 
ambiguities of various terms is difficult to see in approximate calculations. 

As concerns the phenomenological applications, the final state interactions in B hadronic 
decays were investigated the last years by means of dispersion relations with respect to 



the momentum squared (s) of the B meson [18|-||21| (recently, this method was applied 



also to K — > 7T7T decay p^|-|36|). These dispersion relations look more familiar, due to 
their formal resemblance with the case of the pion form factor. However, as seen from Eq. 
(§), the similarity is not complete, due to the presence of the initial strong interactions in 
the weak decay and the appearance of off-shell quantities. The dispersion relations with 
respect to the momentum squared k\ , which seem less intuitive, were applied to B decays 



in Refs. (we mention also one earlier application of this technique for the calculation 
of the nucleon form factor f31j ). In Section 3 we shall use this type of dispersion relations, 
written in the convenient form (fT^), for discussing the effects of the final state interac- 
tions in the B° — > 7r + 7r _ decay. Before making this analysis, we shall first investigate the 
spectral function (Jfsi appearing in this relation. 



3 The Goldberger-Treiman procedure 

Following Ref. [[J, we shall use the parametrization 



A(B° d -> 
T U (B 



i^m 2 Bf+( m l)U\V ud V* b \ e 



1 7 



X 



-17 



7r + 7r~) + -—T C (B -> 7T + 7T-) 



R 



(15) 



obtained by extracting from the amplitudes A u and A c of Eq. (|3|) the "naive" factorized 
amplitude, expressed in terms of the pion decay constant f n and the B — > ir transition 
form factor /+(m£) ( R b = \V ub /(XV cb )\(l - A 2 /2) » 0.377). The spectral function a FS i 



defined in Eq. ( |T0"D can be written in a similar way as 

.Gf 2 



-17 



a 



FSI 



+ 



Rh 



-a 



FSI 



(16) 
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where, recalling that crpsi is given by the first term in the unitarity sum (§), we have 



a FSI ~ J2 S ( k ± + ^2 ~Pn) (7T l^7r+|n)(n|?7^|0), 
n 

j = u,c. (17) 

We assumed here that the source r/£ admits a decomposition in two terms, analogous to 
that of the weak Hamiltonian (fj). This shows also that one can derive separate dispersion 
relations for each of the amplitudes T u and T c . In particular, according to Eq. (|T4|), we 
shall write these relations in the form 



T,=T,, + ^ln 

7T 



9 

ml 



(m B - m T: 



\2 



- 1 



J = u,c, 



where we denoted by T J)0 the analog of the term A in the relation flT4|), divided by the 
constant factorized in ([15]). 

For further applications, it is important to notice that the spectral functions cr FSI and 
cy pgj are real quantities: 

oisj = (^5/)*> J = u,c. (19) 
The proof of these equalities is based on the properties of the matrix elements in Eq. (P~T| ) 



under the PT transformation [30|. More precisely, in the present case, we have 



(it (k 2 ) I r]^+(0)\n, in) = (vr (k 2 )\i] 7T +(0)\n, out)* 

(n,in|4(0)|0) = -(n,out|4(0)|0)*, (20) 

where the minus sign in the second relation is due to the specific spin-parity properties 
of the relevant part of the weak hamiltonian. Using the equalities (|20| ) in Eq. (|TTD, and 
taking into account the equivalence of the complete sets of "in" and "out" intermediate 
states in the unitarity sum, it is easy to prove the relations ( |T9"D (the minus sign in Eq. 
( p0|) is compensated by leaving aside an imaginary constant in the definition (|16|)). 

In approximate calculations, the set of intermediate states in the unitarity sum flTT| ) 
is truncated, which might lead to violations of the reality conditions (0), and to the 
appearance of artificial strong phases in the spectral functions. This fact is important 
in the present case, since the unitarity sum is evaluated at large cm. energy squared, 
s = m 2 B , where many inelastic channels are open. However, with a "good" choice of 
the truncated set one can avoid the appearance of unphysical phases. The idea of Gold- 



berger and Treiman |37| was to take the intermediate states in the symmetric combination 
l/2\n, in) (n, in| + l/2\n, out)(n, out|, 

which represents also a complete set. The remarkable point is that, even when it is 
truncated, this set generates spectral functions which satisfy the reality condition (|T9"D, 
at each step of approximation. The symmetric summation simulates therefore, in a cer- 
tain measure, the effects of inelastic states, without incorporating them explicitly in the 
unitarity sum[]. 

3 An alternative approach to include the effects of the inelastic states in B hadronic decays, based on 
statistical arguments and Regge phenomenology, was proposed in Refs.pq], p7| 
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It is worth mentioning that a complete set written in a symmetric form is quite natural 
in the LSZ method: indeed, when deriving the discontinuity of the amplitude, the initial 
6{xq) function in Eq. (^), whose origin is the reduction of an "out" pion, is actually 
replaced by 8(xq)/2 + 8(—xq)/2 This means that the final two-pion state appears 



in the discontinuity in the symmetric combination l/2|7r + 7r~, out) + l/2|7r + 7r~, in) 4 . It 
is therefore reasonable to take the same symmetric combination also for the intermediate 
states n. 

The Goldberger-Treiman procedure allows us to write the spectral functions (Jpsi 
defined in Eq. (|T7[) as 

4s/ = \ E <K fc i + h - Pn)[M*(n -> 7r + 7r-)3>(S -> n) 

^ n 

+M{n^n + n-)T*(B -> n)], j = u,c, (21) 

where M(ra — > 7r + 7r^) denotes the amplitude of the strong transition from the intermediate 
state n to the final tt + tt~ state, and Tj(B —>■ n) is the specific part of the weak decay 
amplitude of B into the same intermediate state (divided by the constant factorized in 
Eqs. (|15D and flT6|)). 

It is known that the strong dynamics at high energies is dominated by multiparticle 
production. However, as argued in ||27|| , the contribution of the multiparticle intermediate 
states in B decay is suppressed by a flavour mismatch between the weak and the strong 
parts of the process. Therefore, only the states composed of two low mass resonances 



are expected to bring an important contribution to the rescattering. As shown in |27| , 
this picture is consistent with the absence of final state interactions in B decays in the 
heavy mass limit ||, since the production of two resonances is expected to vanish at high 
energies. 

In the two-particle approximation of the unitarity sum, the weak decay amplitudes 
are completely specified by the masses of the particles in the intermediate states, being 
independent on the phase space variables [P^| . Then the phase space integration implicit 
in Eq. (511) can be performed exactly, leading to 



FSI = \ E WtxiPaPa - ^-)T}{B - P a P a ) 



PaPa,X 



+M ,x(PaP a -> 7r + 7r-)T/*(S -> P a P a )} + . . . , (22) 

where M 0i \(P a P a — > 7r + 7r~) denote the S'-wave projection of the strong amplitudes and a 
summation over the helicities A of the intermediate states must be peformed in general. 
By inserting this expression in Eq. ([Tj|), and writing explicitly the real and the imaginary 
part of the logarithm, we obtain the relation 

Tj(B -> 7T + 7T-) = T jfi (B -> 7T + tO 



+ 



E Re[M* >x {P a P a -> n + n-)T 3 x (B -> P a P a )] + ... 



i A similar symmetric combination of "in" and "out" states was obtained in a related context in Ref. 
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i + - In 1 

7T \ 



(m B - 



J 



(23) 



The sum includes the dominant two-particle quasielastic states and resonances P a P a , and 
the dots represent the contribution of the multiparticle states. 

The strong amplitudes entering Eq. ( |23"D are evaluated at the cm. energy squared 
equal to m\ « 25 GeV 2 . For low masses of the intermediate particles we can use the 
generic Regge amplitude p9| 



-7(*) 



t + e 



sin ira(t) 



s \ Q W 

So- 



(24) 



where j(t) is the residue function, r the signature, a(t) = a + a't the linear trajectory, 
and So ~ 1 GeV 2 . Modifications of the standard Regge expression when the particles have 
larger masses are discussed in [27]. The 5-wave projection of the amplitude ([24] ) in the 
spinless case is 



M 0fi (P a P a -> 7T + 7T-) « e 



1 



7(0) 

327ra' rriBqL 



xe 



(a +a't +2aW) i 



-2a'qq'L 



(25) 



where g = l/2^m B — 4m\ an d = 1/2 
and intermediate state, respectively, to 



m\ — 4ml are the c.m.s. momenta of the final 
(q 2 + q' 2 ), L = lnm^/so — in/2 and ^ is a 



signature factor (equal, in particular, to —1 for the pomeron and iy2 for the p trajectory 
25|). In deriving Eq. (f25|) we neglected the t dependence of the ratio j(t)/ sin(7ro;(t)/2) 
for r = 1 trajectories, and of the ratio j(t)/ cos(7ra(t)/2) for r = —1. 



4 Quark-hadron duality 

The relation between the QCD predictions and the hadronic physics is an extremely com- 
plex, still unsolved problem. For testing quark-hadron duality it is in principle necessary 
to perform an analytic continuation from the spacelike region of momenta, where Op- 
erator Product Expansion and perturbative QCD are valid, to the timelike axis, where 
the physical processes are described in terms of hadronic degrees of freedom. This proce- 
dure, based on dispersion relations, was applied to simple objects, like the current-current 
vacuum correlation functions or the electromagnetic form factors. 

The weak hadronic decays are much more complicated, due to the presence of hadrons 
in both initial and final states. Usually, the strong processes at scales larger than nib are 
integrated out, being included in the Wilson coefficients entering the effective hamiltonian 
(0). In perturbative QCD, the decay amplitudes are treated in the heavy mass limit 
using the framework of perturbative factorization for exclusive processes, based on hard 
scattering kernels and light-cone distribution amplitudes, with the heavy mass playing 
the same role as the large momentum transfer. 

In the hadronic picture, we shall consider the dispersion relation (^) (written in 
more detail in (|23"D) where, as discussed in Section 2, the first term T j is given by the 
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contributions which remain after switching off the final state strong interactions among 
the emitted pions. This term should be provided, in principle, by a nonperturbative 
calculation, which excludes in a systematic way the final state interactions. As such a 
calculation is lacking, we resort to a qualitative discussion, based on quark-hadron duality. 

We recall that, according to the discussion below Eq. (0), the terms T j include the 
degenerate terms and the last two dispersive integrals in the relation (|TT|). In all these 
terms, the two final pions appear in different matrix elements, and may be associated 
qualitatively to diagrams with no gluon exchanges between them. In particular, in the 
spirit of the dispersive formalism, the last two integrals in Eq. (|TT| ) may be interpreted as 
"initial state interactions", in a crossed channel. Therefore, it is reasonable to assume that 
a considerable part of T^q consists of the naive factorized amplitude, which is associated 
to processes with no gluon exchanges between the emitted pion tt + and the system (n~B®) 
(except those already included in the Wilson coefficients and the short distance processes 
taking place before hadronization) . As for the last term in the relation (|14D (or (|23|)), 
which describes in the hadronic picture the final state strong interactions, it is dual to the 
topologies (penguin annihilation, exchange diagrams, scattering of the spectator quark, 
vertex corrections to the emission diagrams), involving gluon exchanges between the final 
pions. Of course, the correspondence between various quark diagrams and the terms 
appearing in the hadronic formalism is not simple beyond the lowest orders of perturbation 
theory. In particular, it is impossible to associate in an univoque way the diagrams 
involving many gluons either to the initial, or the final state interactions. Moreover, from 
the point of view of final state interactions, the spectator quark does not play a special role. 
Therefore, the vertex corrections involving the quarks emitted in the weak process, which 
are included in the factorized part of the amplitude in the standard QCD factorization 
approach ||, contribute also to the final state interaction part of the amplitude. 

In the QCD factorization approach ||, the dominant contribution to the decay ampli- 
tude is given by the factorized term, with corrections which are suppressed, in the heavy 
limit nib — » oo, either by powers of a s (mb), or by powers of AQco/^ib- As discussed 
recently, the power suppressed corrections might be enhanced by pure soft effects, such as 
endpoint singularities and higher twist terms in the pion distribution amplitudes, appear- 
ing in annihilation diagrams or in the hard spectator interactions. A recent evaluation of 



these corrections in perturbative QCD factorization was done in but there are still 
differences between the results obtained by different authors. 

In what follows we shall estimate the heavy mass corrections produced by the final 
state interactions, using the hadronic dispersion relation written in the form (f23l). We 
shall consider first the contribution of the elastic and quasielastic rescattering, taking 
as intermediate states P a P a the lowest pseudoscalar mesons 7r + 7r~, 7r°7r°, K°K°, K + K~ 
and 7777. The effect of higher resonances describing inelastic rescattering will be discussed 
below. 

We assume that the amplitudes Tj appearing in Eq. (|3|) can be expanded in the 
heavy quark limit, i.e. for large mg~ m;,,as 

Tj « T jfi + 0(a s (m B )) + 0(A/m B ) + . . . , (26) 

where Tj t o are approximately given by the factorized amplitudes, with short distance 
corrections. These values have small imaginary parts, produced by the complex effec- 
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tive Wilson coefficients, vertex corrections or short distance effects in the penguin and 
annihilation diagrams. 

It is easy to write down the high energy limit of all the quantities entering Eq. (|23|). 
First, from the explicit expression of the logarithm in the r.h.s. of this relation, it follows 
that the real part of FSI amplitude is suppressed by two powers of the heavy mass with 
respect to the imaginary part. The heavy limit behaviour of the Regge amplitude (|25|) 
depends of the specific trajectory. For the pomeron (with ao ~ 1.0 and a' ~ 0.25) we 
obtain the expansion 



M, 



7p(0) 

47T 



7T 



1 



ln^L 

so 



2 i n 2 ™B. 

so 



(27) 



which shows that the dominant contribution in the heavy mass limit is imaginary. At the 

r(P) ~ -0.23 + 0.691 
a 1, we obtain 



physical scale, using 7p(0) ~ 25.6 [^5 



we obtain Mg ~ 

For a physical trajectory, like p, using ao ~ 0.5 and a' 



M, 



(p) 



7,(0) 
WTrrriB 



1 7T 1 



In- 



2ln 2 ^i 

SO 



+ 



:28i 



This amplitude is suppressed by one power of m B compared to the pomeron amplitude 
(f2~7p. At the physical scale, using 7 P (0) ~ 31.4 p5| , we obtain for the p trajectory 
Mq ~ 0.015 + 0.047z. The masses of the particles undergoing the strong scattering appear 
as power suppressed terms in the expansions (|27D and (p8|). 
By inserting the expansions (^ _ 



- (|28|) in Eq. (|23|), we can derive iteratively the 



magnitude of the coefficients of the logarithmic and power corrections in the heavy limit 
expansion ( ^B]) . As an illustration of the method, we take as input of the iterative proce- 
dure the values 



T c ,o 



P< 



Rb 



0.969 - 0.017i 
0.246 + 0.03i , 



(29) 



which are typical for the factorized amplitude with short distance corrections [14 



Using these values as the lowest approximation of the amplitudes Tj in the right hand 
side of the relation (|2"3]), we obtain to first order 



T = T U)1 
Ri, 



P 



0.969 - 0.23i 
0.246 + 0.00012z. 



(30) 



The dominant corrections are given by the pomeron contribution to the elastic channel. 
In particular, the imaginary part of T Uj i is due mainly to the next to leading order (sec- 
ond) term in the expansion ( p7[ ) of the pomeron amplitude, since the dominant term is 
suppressed by the symmetric Goldberger-Treiman summation in Eq. (^3|). 

The values fl3"UD have uncertainties due to the higher order terms in the heavy mass 
expansion. First we notice that the contribution of the pseudoscalar mesons 7r°7r°, K°K°, 
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K + K~ and 7777, responsible for the quasielastic rescattering in the unitarity sum of Eq. 

is negligible. Indeed, these states are produced by non dominant decay diagrams, 
and the corresponding Regge amplitudes are described by physical trajectories, which 
are suppressed with respect to the pomeron. The estimates made in Ref. [|25|| , based on 
SU(3) flavour symmetry, show that the effect of these channels on the spectral functions 
is not larger than several percents. 

The inelastic states might bring an important contribution in the unitarity sum if 
they are produced by dominant decay diagrams, or the corresponding CKM coefficients 
are large. The first states which can contribute are the lowest vector resonances, p + p~, 
p°p°, K*K*, uouo and 00, which describe intermediate states with four or six pions. The 
quark diagrams for B decays into VV pairs are similar to those of the corresponding 
pseudoscalar mesons. Therefore, only the pair p + p~ is produced by a dominant tree 

M 



give for the helicity amplitudes of the 
up to a factor of 2, with the amplitude 



diagram. Estimates based on factorization 
weak decay — > p + p~ values comparable 

of the decay B% — > 7r + 7r~. As concerns the strong amplitude M(p 
described in the Regge model by the exchange of uo and A 2 , with trajectories a ~ 0.5 + 1 
(almost degenerate with that of p), and the 7r exchange, whose trajectory a w ~ 0. + 0.9t 
is non-dominant [|39|. As discussed in Ref. Ell, the kinematic factors of the t-channel 



p — > 7r + 7r ), it is 



J9| . As discussed in Ref. 
helicity amplitudes suppress in the present case the contribution of the natural parity 
exchanges u and A 2 , so that the scattering amplitude of p + p~ — > 7r + 7r~ is described at 
small t by tt exchange. This implies a considerable suppression, verified indirectly in the 



related reactions of vector meson production tc N — > p N [[12|, ^3[. Therefore, the vector 
meson resonances seem to bring a negligible contribution to the inelastic rescattering 
in B® — > 7r + 7r~ decay. Using this spin suppression argument, we expect that the same 
conclusion applies also to other vector mesons of higher mass. 

The role of the intermediate states with charm, like the pair DD, has been considered 
by several authors and is still controversial. Since a large part of the inelastic irn scat- 
tering at ^/s ~ tub goes into multiparticle states composed of noncharmed mesons, the 
contribution of these states was assumed in Ref. [f?7| to be negligible. The annihilation 
of the cc pair in the penguin diagrams proceeds therefore through a short distance inter- 
action, which can be computed perturbatively. We mention however that the possibility 
of a long distance contribution of the "charming penguins" was also considered by some 
authors | 44fl , possibly through the intermediate state DD 



23L 24 



A detailed estimate is 

difficult since the validity of the Regge model is questionable at yfs ~ m B for large masses 
(mu = 1.68 GeV) of the particles undergoing the strong rescattering [p7[| . We mention 
that the contribution of high mass resonances, such as DD, is actually suppressed by the 
phase space integral in Eq. ([H]), evaluated at fixed s = m B . In the present approach, 
we assume also that their global effect is taken into account in a certain measure by 
the Goldberger-Treiman symmetric summation. With this assumption, we do not expect 
other important power corrections in the spectral functions. 
From the values given in Eq. ([SO) we obtain 



0.246 , 5 = Arg 



14° 



(31) 



while the input values ( p9|) correspond to |P /^o| = 0.25 and Arg[P /To] = 8°. Thus, in 
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the present approach the final state interactions do not modify considerably the modulus 
and the phase of the ratio P/T. We recall that the lowest order values Tj t0 depend on 
the renormalization scale \x in the Wilson coefficients. It is expected that this dependence 
will be diminished by the inclusion of other terms in the l.h.s. of Eq. (E£3). 

5 Conclusions 

In the present work we investigated the effects of the final state interactions to B® — > 
7r + 7r~ decay in a formalism based on hadronic unitarity and dispersion relations. We 
discussed the heuristic derivation of the dispersion relations with respect to the momenta 
of the external particles, by applying the LSZ procedure to the S-matrix element of the 
weak decay. The ambiguities which affect in general the off-shell extrapolation of the 
amplitudes, appear in our formalism in the so-called "degenerate terms", produced by 
the equal-time commutators in the LSZ formalism, and in the off-shell quantities entering 
the spectral functions. 

The evaluation of the dispersion relations in external momenta is in general very com- 
plicated. However, we noticed that the dispersion relation with respect to the momentum 
squared of one final pion can be written in the convenient form (|T4]) , where the contri- 
bution of the final state interactions is separated from other terms, and involves only 
on-shell quantities. We used this relation, written in more detail in Eq. (^), as an itera- 
tive scheme for determining the corrections to the factorized amplitude, generated by the 
final state interactions in the heavy mass limit. A nontrivial prediction of the formalism 
is that the real part of the FSI contribution is suppressed by two powers of the heavy 
mass, compared to the imaginary part. Using for illustration a numerical input suggested 
by QCD factorization, we noticed the dominant effect of the next to leading logarithmic 
term of the pomeron contribution. Other sources of large power corrections to the factor- 
ized amplitude are not found. We discussed the contribution of the lowest pseudoscalar 
mesons and vector meson resonances, and assumed that the effects of higher resonances 
and multiparticle states are qualitatively taken into account by the Goldberger-Treiman 
method of calculating the spectral functions. In particular, the results of a numerical test 
indicate that the phase and the modulus of the ratio P/T are not drastically modified by 
the final state interactions. Using improved results of QCD calculations, it will be possible 
to test the dispersion relations conjectured in the present work and, more generally, the 
validity of quark-hadron duality. 
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